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Abstract 



We perform a Monte Carlo simulation calculation of the critical coupling 
constant for the continuum jt 



b\ theory. 



The critical coupling constant we 



obtain is 



A 

J?. 



crit 



10.24(3). 
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1 Introduction 



The 02 fi e ld theory, specified by the (Euclidean) Lagrangian 



(i) 



has solutions in a symmetric phase in which the discrete symmetry of the Lagrangian, 
(j) — > —(f), is manifest, i.e. < >= and there is no trilinear coupling. It also has 
solutions in a broken symmetry phase with < (f> >^ and induced trilinear couplings 
proportional to < >. 

There exist both elegant heueristic |l| and rigorous || mathematical proofs of the 
existence of this phase structure, but there is no rigorous result for the critical value 
of some coupling constant which separates these two phases. There do exist numerous 
approximate calculations, and in the case of the lattice model (lattice spacing a > 0), 
the critical line in the //jj, A plane which separates the two phases is known to some 
numerical accuracy by Monte Carlo simulation. 

For the Euclidean Quantum Field Theory (EQFT), which is the continuum limit 
(a — > 0) of the lattice model, the first step is to specify what finite dimensionless 
coupling constant is to be used to parametrize the solutions and for which we are to 
determine the critical value. 

For the lattice model, these considerations are straightforward. There are two 
parameters, /Xq, A in Ce, and there is the lattice spacing a. In d = 2 both ^ and A 
have dimension mass squared. (We assume the infinite volume limit, L — > oo). So 
there are two independent dimensionless parameters which may be taken to be the 
two Lagrangian parameters measured in units of inverse lattice spacing squared, 



In this parametrization,the phase diagram of the lattice model consists of the critical 
line in the , Xi plane, determined by Monte Carlo simulation (Fig. 1). 

The situation for the EQFT is more complicated. There is an infinite mass renor- 
malization. The bare mass parameter has to be tuned to infinity as the continuum 
limit is taken, ~ fi 2 In-, where jj? is some finite renormalized mass squared. 
Thus (Xq cannot be used as a parameter of the continuum solution. In d = 2, the 
field strength and coupling constant renormalizations {Z§, Z\) are finite and can be 
disregarded in the study of the phase structure of the theory. Furthermore, since the 
dimensionful coupling constant A is independent of a, and fi^ diverges only logarith- 
mically with a, both A^, /j,q go to zero in the continuum limit, a — > 0. That is, the 



\ \ 2 2 2 2 

Xi = Xa , fj. = fj, a . 
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Figure 1: Phase Transition Line in the A;,^ Plane 

EQFT limit is the single point at the origin of the fJ% g ,\e plane in Fig. [l|. Taking 
the limit a — > reduces the number of of independent dimensionless parameters from 
two to one. 

The required mass renormalization can be written as a simple reparametrization 
of C E (TJ). Let 

al = M 2 - V- (3) 

Then 

£ E = i(V0) 2 + ^ 2 2 + ^ 4 -^ 2 2 (4) 

There is still a substantial freedom of choice in the definition of the finite renormalized 
mass parameter, /A The ultraviolet, ln^, dependence of /i 2 , is moved entirely to the 
counter term 5fi 2 ; but the separation of the finite part of //q into /i 2 and 5fi 2 is only 
determined when a renormalization condition is specified. The dimensionless effective 
coupling constant = \ then manifestly depends on the choice of renormalization 
condition which fixes the finite part of Sfj 2 . For example, we could take fx 2 = m 2 , the 
pole mass, by choice of renormalization condition 

= G-\p 2 = -ml) (5) 

and dimensionless coupling 

9 = A (6) 



m 2 



A closely related choice, convenient for lattice Monte Carlo measurement, is to take 
/i 2 = m' 2 , defined by the renormalization condition 

m ' 2 = V = 0) (7) 
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and 

9' = A- (8) 

In fact, neither of these choices provide a dimensionless coupling constant whose value 
distinguishes between the two phases of the theory - because the renormalization con- 
ditions themselves do not distinguish between the two phases. Either renormalization 
condition, (|5]) or (|7]), can be implemented perturbatively in either the symmetric phase 
or the broken symmetry phase. The phase has to be specified by ansatz, < <fi >= 
or < <fi >^ 0, so that one can perturb about the correct stable vacuum. Then g or g' 
can take on arbitrarily small values in either phase. 

A dimensionless coupling constant whose critical value separates the two phases 
is provided by choosing the mass renormalization to be equivalent to normal ordering 
the interaction in the interaction picture in the symmetric phase. From (fj]),([|) 

G- V) = P 2 + Mo + S (P 2 ) = V 2 + ^ + W) (9) 

and for /i 2 > 

S(p 2 ) = 3 A - b\x 2 + two-loop. (10) 
A^2 in the continuum limit is the ultraviolet divergent Feynman integral 

A„2 = / — 11 



(2^)V + /i 2 ' 



On the lattice, 



^ N N | 



A " 2 £ A + 4 sin 2 2f + 4 sin 2 ^ (12) 

The 'leaf diagram (Fig. 2) which gives the integral A^2 is the only divergent Feynman 
diagram of the theory in d = 2. 



Figure 2: Leaf Diagram 
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Thus the renormalization condition 



5/i 2 = 3 A A^2 (13) 

removes all ultraviolet divergence from the perturbation series based on the renor- 
malized parametrization given by (|j) and ( jT3[) . 

= ^(V0) 2 + ^V + ^0 4 V (14) 

In the last normal ordered form, we have dropped a constant piece. The dimensionless 
coupling constant suggested by (|TJj) we denote simply by /. 

f=~ 2 (15) 

Since the normal order prescription in (FJ) is relative to the vacuum of the symmetric 
phase theory, we may investigate the possibility of a critical value of the coupling 
constant /. Using ([15]), the first line of ( |H| ) may be rewritten as 

£e = i(V0) 2 + ^ 2 (1 - SfA^ 2 + (16) 

On the lattice (fixed a > 0), A„2 is finite; and we can argue that for small enough 
/, the exact effective potential is well-approximated by the classical effective poten- 
tial with its single minimum at (p c i = 0. For large /, the coefficient of 2 in ( |3~6| ) 
becomes negative, suggesting a transition to the broken symmetry phase. However 
the argument falls short at this point because for strong coupling one can not argue 
that the effective potential is well approximated by its tree level form. The argument 
was completed by Chang H by constructing a duality transformation from the strong 



coupling regime of (|14] ) to a weakly coupled theory normal ordered with respect to 
the vacuum of the broken symmetry phase. 

There are several attempts in the literature (see Table 2) to compute the critical 
value, f c , by various analytic approximations, with a rather large spread of answers. 
In this paper we report an accurate numerical value by Monte Carlo simulation. The 
first step is to obtain the critical line in the Hq , plane. This determines ^Q,{\t) C rit 
(see Fig.l). Recall that these values are infinite volume extrapolations of finite volume 
Monte Carlo data. Then, combining (|3|),(|T3|) we obtain 



= ~ 3 X t A » 2 (I 7 ) 
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In the infinite volume limit A^i fll^) has the integral representation 

/•CO 

An2 = / dt exp(-/^ t)(exp(-2t)J (2t)) 2 
jo 



For any point away from the origin, (|17| ) , ([Tq) can be solved numerically to determine 
crit' This is then extrapolated into the origin to determine 

f c = lim -|| cri t (19) 

2 Simulations 

The Monte Carlo simulation is based on the lattice action which regularizes the 
continuum theory (|1|). 

^ = EUE(^+^)-^)) 2 +\vlrtn) 2 + jV(nA (20) 

Periodic boundary conditions were imposed on N x iV square lattices with iV = 
32, 64, 128, 256, and 512. To reduce critical slowing down, our updating algorithm 
consisted of a standard Metropolis update (i.e. with a symmetric transition matrix) 
alternating with a cluster algorithm updating the embedded Ising model. The pro- 
cedure is similar to that of Brower and Tamayo ||, but we substitute a Wolff- type 
single-cluster algorithm for the Swendsen-Wang multiple-cluster algorithm used there. 

Measurement of lattice quantities was performed every ten Metropolis+cluster 
cycles. Each data collection run consisted of 10 4 to 10 5 measurements, after an initial 
thermalization of at least 10 4 cycles. To assess the effective number of statistically 
independent measurements, the integrated autocorrelation time (tint) was calculated 
for each run, 

1 ^ a(i) , , 

where s(z) is the autocorrelation separated by i measurements and M is some num- 
ber of measurement such that s(M) is essentially noise. The largest tint measured 
in our simulations was 8 measurements (80 update cycles) and was typically much 
smaller. Thus in every case the thermalization time exceeded IOOtint, so we expect 
that our lattices were well-thermalized before we began collecting data. We also tried 
to measure the exponential autocorrelation time texp from the first few (time) auto- 
correlation functions, but these didn't appear to fall off as a single exponential. As 
expected, however, the measured values were slightly smaller than the corresponding 
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tint for that run. As a result of the small tint and large number of measurements, 
statistical errors are typically quite small, generally smaller than the systematic errors 
in the determination of the critical line in the A^, y? h plane for finite volumes and in 
the extrapolation to the infinite- volume limit. 

For each size lattice, we did looked at Xe = 1.0, 0.7, 0.5, 0.25, 0.1 and 0.05, and for 
each Xi scanned in (jl'q , starting in the symmetric phase and ending in the broken 
symmetry phase. We used two diagnostics to determine the critical value of /Xq . The 
first was the value of /Xq which produced the maximum of the variance of the action 
(specific heat), which should diverge as Ch ~ In I/Xq — /Xq | as lattice size L — > oo [Q. 
The second was based on the shape of the histogram of the distribution of < </? >. 
In the symmetric phase with L > ^, the spatial correlation length, the probability 
distribution of < (p > should be a single peak centered about zero, while in the 
broken-symmetry phase with L ^> £ the distribution should consist of two identical 
peaks at equal distance from < <p >= 0. For a fixed value of /Xq , histograms of 
< ip > will approach one of these distributions for L sufficiently large. To quantify 
the bimodality of the distribution, we define h to be the ratio of the number in the 
central bin to the largest number in any outlying bin. For a histogram which is a 
single peak centered about zero, h > 1. For a histogram which is two peaked, around 
+ < \<p\ > and — < \<p\ >, h « 1. The diagnostic for a given /Xq is the behavior of h 
as the lattice size is increased. For /Xq 's which lead to symmetric phase in the infinite 
volume limit, h increases with L until it exceeds one. For /x^'s which lead to broken 
symmetry phase in the infinite volume limit, h is rapidly decreasing with increasing 
L. For a narrow range of /Xq around the infinite volume critical value, this behavior 
may not stand out until one gets to quite large lattices. In practice, on the lattices 
we used, from these observations we obtained upper and lower bounds on the infinite 
volume critical value of /Xq which are closer togther than the estimated statistical 
and systematic error in the value obtained from extrapolation of the position of the 
peak of the variance of the action. (In all cases, within that estimated error the 
extrapolated value is consistent with the bounds). 

3 Analysis 

In the second column of Table 1 we give our estimates of the critical value of /Xq for 
each of the values of Xe listed above, extrapolated to the infinite volume limit. (They 
are the input for Fig. 1). They may be compared with the values found by Toral and 
Chakrabarti 0] seven years ago, extrapolating from much smaller lattices. Within 
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A 


f'OC 




A 


1.0 


-1.270(3) 


0.0980(8) 


10.204(80) 


0.7 


-0.9516(8) 


0.06844(23) 


10.228(33) 


0.5 


-0.7210(10) 


0.0489(3) 


10.225(63) 


0.25 


-0.4035(5) 


0.0242(2) 


10.33(8) 


0.10 


-0.1878(5) 


0.009615(140) 


10.40(15) 


0.05 


-0.0998(3) 


0.00492(9) 


10.16(19) 



Table 1: Determination of the Phase Transition Line for Different A; 




Figure 2: Plot of Phase Transition Line in the /x 2 , \i Plane 

the estimated errors they are consistent, except for the smallest A^, which is closest 
to the continuum limit and requires larger lattices to accomodate both large L and 
small a. In the third column of Table 1 we give the corresponding critical values of /j 2 
as determined from fll^D , (|18D, and in the fourth column we give the corresponding 
values of 

In order to get a feel for the systematic errors we have done the extrapolation to 
the continuum limit in a number of different ways. We have plotted [i 2 vsA and taken 
the inverse of the slope at the origin to determine the critical value of /. This is 
shown in Fig. 3. Note that the data fall very well on a straight line (x 2 = 0.65) which 
passes throught the origin within the estimated error. As a small variation on this we 
have redone the fit forcing the line to go exactly through the origin. The two values 
of f c determined in this manner are 10.23(3) and 10.24(2) respectively. Alternately, 
we have extrapolated the values of / = to the continuum limit (A; ai — > 0) as shown 
in Fig. 4. We assume that we approach the continuum limit along a 'line of constant 
physics', i.e. once we are close enough to critical (£/ at ^> 1), dimensionless physical 
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Figure 3: Extrapolation to determine the critical value of / 

quantities should not change their values as the lattice spacing is reduced. Thus we 
fit the data in Fig.|] with a horizontal line. The \ 2 = 0-59 is consistent with this 
assumption and the result for f c is 10.24(3). 

Although the errors originated as statistical errors in simulations on finite lattices 
the subsequent extrapolations have introduced systematic errors larger than the sta- 
tistical errors. Thus the assignment of the final result and errors are just based on 
the consistency of the above numbers. We conclude that the critical value of is 
10.24(3). 



4 Analytic Approximations 

The above result may be compared with various approximate analytic calculations. 
The simplest approach would be to consider the one-loop effective potential. In d = 2 

tilllS IS 

^/ = y 2 + I 4 -^^(ln— ^ 1) (22) 

With the assignment k 2 = fi 2 this V e ff gives a first-order phase transition for f c = 6.6. 

A compendium of nonperturbative analytic approximate calculations has been 
compiled in Ref. ||. We include these results for purposes of comparison with the 
numerical MC simulation result. 

There is also the issue of the order of the phase transition. According to the 
Simon-Griffiths theorem |TTJ, the phase transition is second order. In the analytic 



approximations the order of the phase transition is determined; the one-loop effective 
potential and the Gaussian effective potential predict a first-order phase transition, 
while the other correctly predict a second-order phase transition. In our numerical 
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Approximation 


Result 


Reference 


Non-Gaussian Variational 


P. 88 
O.oo 




Discretized Light-Front 


7.316,5.500 


EM 


Coupled Cluster Expansion 


3.80 < f c < 8.60 


E 


Connected Green Function 


9.784 


E 


Gaussian Effective Potential 


10.211 


1 




10.272 


i 



Table 2: Analytic Approximations of the Criticaf Value of / 



calculation of the of the critical couping constant we have made no effort to distinguish 
between weakly first-order and second-order phase transitions. 

5 Conclusions 

We have calculated an accurate numerical value of the critical coupling constant us- 
ing Monte Carlo simulation. With this we can evaluate the accuracy of analytic 
approximation methods. It is interesting to observe that the Gaussian effective po- 
tential result for the critical coupling is consistent with the accurate numerical result, 
although it gives incorrectly the order of the phase transition. 
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